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ABSTRACT 

We present the spectrum of compact object masses: neutron stars and black holes that originate from 
single stars in different environments. In particular, we calculate the dependence of maximum black 
hole mass on mctallicity and on some specific wind mass loss rates (e.g., Hurley et al. and Vink et 
al.). Our calculations show that the highest mass black holes observed in the Galaxy Afbh ~ 15 M© in 
the high metallicity environment {Z = Zq = 0.02) can be explained with stellar models and the wind 
mass loss rates adopted here. To reach this result we had to set Luminous Blue Variable mass loss rates 
at the level of ~ 10~^ Mq yr~^ and to employ metallicity dependent Wolf-Rayet winds. With such 
winds, calibrated on Galactic black hole mass measurements, the maximum black hole mass obtained 
for moderate metahicity {Z = 0.3 Zq = 0.006) is Mbh,max = 30 Mq. This is a rather striking finding 
as the mass of the most massive known stellar black hole is Mbh = 23 — 34 M© and, in fact, it is 
located in a small star forming galaxy with moderate metallicity. We find that in the very low (globular 
cluster-like) metallicity environment the maximum black hole mass can be as high as Mbh, max = 80 Mq 
[Z = 0.01^0 — 0.0002). It is interesting to note that X-ray luminosity from Eddington limited accretion 
onto an 80 M© black hole is of the order of ^ IQ^^ erg s~^ and is comparable to luminosities of some known 
ULXs. We emphasize that our results were obtained for single stars only and that binary interactions 
may alter these maximum black hole masses (e.g., accretion from a close companion). This is strictly 
a proof-of-principlc study which demonstrates that stellar models can naturally explain even the most 
massive known stellar black holes. 

Subject headings: binaries: close — black hole physics — gravitational waves — stars: evolution — stars: 
neutron 



1. INTRODUCTION 

Measuring the masses of celestial objects is one of the 
principal challenges in astrophysics. The primary method 
used to accomplish this task involves considering orbital 
motion. This standard method has been applied to a num- 
ber of X-ray binaries with the result showing that the mass 
function exceeds the maximal mass of a neutron star (NS), 
i.e., 2 — 3 M0, which points to the fact that there arc sys- 
tems containing black holes (BH). There are also indirect 
ways of measuring black hole masses. The measurement 
of X-ray luminosity leads to a lower limit on the mass of 
the accreting body arising from the Eddington limit. This 
method, which was applied in the case of ultra luminous 
X-ray sources (ULX) has hinted that these systems are 
BHs with masses in excess of 100 Mq (Miller et al. 2004). 
This claim is subject to a number of assumptions, such as 
the isotropy of radiation (e.g.. King 2009). Additionally, 
it was implied that BHs with mass smaller than 80 M0 
can explain all, but the most luminous, ULXs (for most 
recent studies see Zampieri & Roberts 2009; Kajava & 
Poutanen 2009; Gladstone, Roberts & Done 2009; Mapelh 



et al. 2009). Nevertheless, there is a growing sample of 
stellar mass black holes with masses confirmed by observa- 
tions of orbital binary motion, for a discussion see Pakull 
& Mirioni (2003). Moreover, there are some solitary BH 
candidates which have been detected through microlens- 
ing experiments, and there are good prospects for observ- 
ing more of such objects. A review of the observations 
of galactic BH candidates (Orosz 2003; Casares 2007; Zi- 
olkowski 2008) shows that the masses of stellar BHs range 
from a few solar masses to about 14±4 M0 (Grcincr, Cuby 
& McCaughrean 2001), in the case of GRS 1915-105, and 
16 ±5 for Cyg X-1 (Gies & Bolton 1986). However, in case 
of Cyg X-1 the estimates by different authors vary from 
10 M0 (Hcrrero ct al. 1995) to 20 ± 5 M0 (Ziolkowski 
2005). 

The development of experimental techniques in recent 
years has allowed for the investigation of X-ray binaries 
in neighboring galaxies. This has led to the discovery of 
the largest stellar mass black hole in the binary ICIO X-1. 
The mass of the BH was found in the range of 23 — 34 Mq 
(Prestwich et al. 2007; Silverman & Filippenko 2008). It 
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is interesting to note that the star formation rate in dwarf 
galaxy ICIO is very high, and its metahicity is low ~ O.SZ© 
(Massey et al. 2007). 

The spectrum of black hole masses is extremely interest- 
ing from the point of view of gravitational wave astronomy. 
The interfcrometric observatories like LIGO (Abramovici 
et al. 1992) and VIRGO (Bradaschia ct al. 1990) can de- 
tect coalescences of 10 Mq BH binaries at a distance of 
« 150Mpc, and the sensitivity of the detectors will soon 
increase. Thus, the theoretical understanding of the dis- 
tribution of masses of black holes formed by stars may 
soon be measurable with gravitational wave observations, 
and vice versa; knowledge of the BH mass spectrum may 
help to identify the parameter space which favors a high 
probability of source detection. 

Formation of BHs in the course of stellar evolution is 
connected with the end of the nuclear burning phase in 
a massive star. For the lower mass end of BH formation, 
a meta-stablc protoneutron star may be formed, and the 
black hole appears after accretion of the part of the stel- 
lar envelope that could not be expelled in the supernova 
explosion. In the case of the most massive stars the BH 
is formed through accretion of the entire stellar material 
(direct collapse or failed supernova). Thus the mass of the 
BH is determined mainly by the mass of the star at the 
moment the collapse takes place, as well as the explosion 
energy (e.g.. Fryer 1999). The presupernova mass is set 
predominantly by the amount of mass loss during stellar 
evolution. 

The mass loss for massive stars is due to stellar winds 
(single stars; e.g., Vink 2008) and additionally close inter- 
actions for stars that are found in multiple systems (e.g., 
binary stars; Hurley, Tout & Pols 2002). Knowledge of 
the mass-loss rates for BH progenitors such as Wolf-Rayet 
(WR) stars - massive stars found near the main-sequence 
(MS) losing mass at high rates and showing weak, or no, 
hydrogen lines in their spectra - and Luminous Blue Vari- 
ables (LBVs) - extremely massive post-MS objects in a 
stage of evolution prior to becoming WR stars - is there- 
fore important for understanding BH masses. A subset of 
WR stars (those with no hydrogen lines) are the naked he- 
lium stars which can also be formed from less massive stars 
that lose their hydrogen envelopes on the giant branch or 
beyond. 

In this paper we review the recent results on modeling 
the stellar winds with a special emphasis on the metahic- 
ity dependence of these winds. This is a proof of principle 
study as we consider only the simplest case: that of single 
stellar evolution. We neglect the effects of binary interac- 
tions (e.g., mass loss/gain due to a close companion) and 
calculate the BH mass spectrum for single stars only. We 
combine the wind mass loss rates with the stellar evolution 
models and investigate the initial-remnant mass relation of 
stars for different metallicities. In particular, we calculate 
the maximum black hole mass for various models. The 
description of the wind mass loss rates and evolutionary 
model is presented in section §2, §3 contains the results 
of our calculations, and we summarize our findings in § 4. 

2. MODEL 

For our study we employ the single star evolutionary 
formulae of Hurley, Pols & Tout (2000) that are used 



within the StarTrack population synthesis code (Belczyn- 
ski, Kalogcra & Bulik 2002; Bclczynski et al. 2008). Up- 
dates to the Hurley et al. (2000) formulae in StarTrack 
include an improved prescription for final remnant masses 
(see below) and the modeling of electron-capture super- 
novae where electrons arc captured onto Mg atoms in an 
0-Ne-Mg stellar core leading to collapse to a NS (Podsi- 
adlowski et al. 2004), thus extending the range of stellar 
masses that produce NSs (at the low-mass end). We also 
note that Hurley et al. (2000) fitted their stellar evolution 
formulae to detailed models of stellar masses of 100 Mq or 
less and we extend this to include stars up to 150 Mq - the 
predicted maximum stellar mass for star formation under 
usual conditions (Weidner & Kroupa 2004). However, the 
formulae are well-behaved within this extrapolation. Stars 
more massive than 150 A/q which are now being consid- 
ered as possible progenitors of intermediate-mass BHs in 
dense star clusters are not considered (see Glebbcek ct al. 
2009 for an overview). 

The employed single star evolutionary formulae were ob- 
tained for stellar models without mass loss. The formulae 
include the effects of mass gain (and potential rejuvena- 
tion) and mass loss on a star. Change of star mass may 
lead (e.g. for main sequence stars) to change of star central 
temperature and pressure, that affects the rate of nuclear 
reactions and a star lifetime. This leads also to the change 
of external star properties like its luminosity and radius. 
There is a feedback running from adopted wind mass loss 
rate (that depends on star properties; mass, luminosity 
and radius) to the star (mass loss affects star mass and 
thus its radius and luminosity) and then back to the wind 
mass loss rate. This setup allows for employment of differ- 
ent wind mass loss rate prescriptions with the same set of 
underlying stellar models. The scheme that is used in our 
calculations is described in detail by Tout et al. (1997) 
and Hurley ct al. (2002). This scheme is only an ap- 
proximation in the mass loss/gain treatment, and detailed 
stellar evolutionary calculations (e.g., Timmes, Woosley & 
Weaver 1996; Limongi & Chicffi 2006) with the specific set 
of winds could be, in principle, used to obtain the more 
accurate results. However, it is noted that both stellar 
models (e.g., mixing or reaction rates; Cassisi 2009) and 
wind mass loss rates (e.g., clumping or LBV phase; Vink 
2008) are burdened with a number of uncertainties, ren- 
dering any estimate of a black hole mass a subject to large 
systematic errors. 

For stellar winds we use both the prescriptions given 
originally by Hurley et al. (2000: sec §2.1) and a new set 
of stellar winds (see §2.2). In the descriptions of these 
winds the following symbols are used: L {[Lq]), R ([i?©]), 
M {[Mq]), and T {[K]) for stellar luminosity, radius, mass 
and effective temperature, respectively, as well as Z for 
mctallicity with the solar value being Z = Zq = 0.02. The 
wind mass loss rates are denoted as dM/dt {[ M© yr^-'^]). 



2.1. Hurley et al. winds: previous reference model 

Here, as it is important for the presentation of the re- 
sults and the discussion, we reiterate the wind mass loss 
prescriptions of the original source (Hurley et al. 2000) 
that have been used in StarTrack over the last several 
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years. The wind prescriptions were adopted as follows: 

(dAf/dt)R, = 2 X 10-'^^ Mq yr-i (1) 

for stars on the Giant Branch and beyond (Kudritzki & 
Reimers 1978; Iben k Renzini 1983); 

\og{dM/dt)yw = -11.4 + 0.0125[Po-100max(A/-2.5,0)] 

(2) 

for stars on the Asymptotic Giant Branch (Vassilidis & 
Wood 1993), with a maximum value of ((iA//(it)vw = 
1.36 X 10~^L Mq yr~^ and pulsation (Mira) period for 
these stars being log(Po) = min(3.3, —2.07 — 0.9 log M + 
1.94 log i?); 

(dA/M)Nj =9.6xlO-i^i?°«iLi-24M°i6 (^)°' ^0 y^'' 

(3) 

for luminous/massive {L > 4OOOL0) stars (Nieuwenhui- 
jzen, H., & de Jager, C. 1990; Kudritzki et al. 1989); 

{dM/dt)^n, = lO-i^Li-S Mq yr-i (4) 

for WR stars (Hamann & Koesterke 1998), it is noted 
that WR-like winds are used also for stars with small H- 
rich mass envelope ( WR star shines through) , although in 
such a case the above formula is modified by the factor of 
1—fJ, where ^ stands for (in first order approximation) frac- 
tional envelope mass (for details see equation 97; Hurley 
et al. 2000); 

(dA/M)LBVi = 0.1(10-^i?LO-^-l)3 (^-^^j^ - 1) Mq yr-i 

(5) 

for Luminous Blue Variables (i > 6 x 10^ and 
10-^RL°-^ > 1.0); Humphreys & Davidson 1994). The 
LBV mass loss rate, if applicable, is added on top of 
stellar winds for H-rich stars, which in turn is calculated 
from {dM/dt) = Ynsx[{dM / dt\, {dM/dt)yw, {dM/dt)^,,]. 
For WR stars the wind is obtained from (dM/dt) — 
max[(dAf/dt)R, (dAf/di)wRi]- 

2.2. Vink et al. winds: new reference model 

For hot massive H-rich stars (B/0 spectral type) follow- 
ing Vink, de Koter & Lamers (2001) we apply: 

\og{dM/dt)oB = -6.688 + 2.210 log(L/105) 

-1.339 log(Af/30) - 1.601 log(y/2.0) 
4-0.85 log(Z/ZQ) + 1.071og(r/20000) 

(6) 

with the ratio of wind velocity at infinity to escape velocity 
V = Voc/vesc = 1-3 for stars with 12500 < T < 22500K; 
and 

\og{dM/dt)oB = -6.697+ 2.1941og(L/105) 

-1.313 log(A//30) - 1.226 log(l//2.0) 
-h0.851og(Z/ZQ) + 0.933 log(T/40000) 
-10.92[log(T/40000)]2 

(7) 

with V = 2.6 for stars with 27500 <T < 50000K. Around 
T ^ 25000K there is a bi-stability jump that leads to 
rapid wind increase. The jump is due to recombination 



of the Fe IV to the Fe HI ion which is a more effective line 
(wind) driver. In the transition zone, we apply eq. [B] for 
22500 <T < 25000K and eq. [7]for 25000 <T < 27500K. 
We note that the above winds apply to stars that are more 
massive than A/zams 3 — 3.5 Mq (spectral type earlier 
than B8-7) as these stars have T > 12500-/^ at the Zero 
Age Main Sequence. 

For Luminous Blue Variable stars that are evolved (be- 
yond the main sequence), H-rich and extremely luminous 
(L > 6 X 10^ and lO'^RL^-^) > 1.0; Humphreys & David- 
son 1994) we adopt 

{dM/dt)i^BV2 = ./ibv X 10-4 Mq yr-1 (8) 

with the standard choice for a calibration factor /ibv = 1.5 
(see §3.3 for justification). This formula is assumed to 
give the full amount of LBV mass loss (unlike in the pre- 
vious prescription in which LBV winds were added on top 
of the underlying massive star winds: § 2.1). We note that 
this formula accounts both for LBV stellar wind mass loss 
as well as possible LBV shell ejections. Furthermore, it is 
assumed that LBV mass loss is independent of metallicity 
(see §3.3 for further discussion). 
For Wolf-Rayet (naked helium) stars we adopt 

{dM/dt)wR2 = 10-"Li-5 l_\ yr-i (9) 

which is a combination of the Hamann & Koesterke (1998) 
wind rate estimate that takes into account WR wind 
clumping (reduced winds), and Vink & de Koter (2005) 
wind Z-dcpendcncc who estimated m = 0.86 for WR stars. 

For H-rich low mass stars, for which the above prescrip- 
tions do not apply, we use Hurley et al. (2000) winds 
(§2.1). 

2.3. Compact object masses: maximum BH mass 

To calculate the mass and type of the compact object 
remnant we use Hurley et al. (2000) single star evolution- 
ary formulae with the updated winds (see § 2.2) to obtain 
the presupernova star mass and structure (i.e., star and 
CO core mass). For a given CO core mass we estimate the 
final FeNi core mass using evolutionary models of Timmes, 
Woosley & Weaver (1996). We use hydrodynamical calcu- 
lations of supernovae explosions to estimate the amount of 
fall back that may occur during core collapse of the most 
massive stars (e.g.. Fryer 1999; Fryer & Kalogera 2001). 
Finally, to change baryonic mass to gravitational mass of 
a remnant we use a prescription proposed by Lattimer & 
Yahil (1989) for neutron stars, and for black holes we as- 
sume that the gravitational mass is 90% of the baryonic 
mass. 

Our prescriptions for the iron core mass, fallback mass 
and final remnant mass are given in full detail in Belczyn- 
ski et al. (2008, see their §2.3.1). For CO core masses 
below 5 Mq, the remnant is set to the iron core mass. For 
CO core masses above 7.6 Mq, we assume the star col- 
lapses directly to a black hole and its final mass is equal 
to the star's total mass at collapse. In between these two 
extremes, we use a linear fit between no fallback onto the 
iron core and complete fallback. This fallback prescription 
is based on the explosion models from Fryer (1999) and is 
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in agreement with Fryer & Kalogcra (2001). We can com- 
pare this prescription directly to the latest analytic study 
of remnant masses by Fryer et al. (in preparation) using 
CO cores from the Heger et al. (2003) solar-metallicity 
progenitors. The maximum remnant masses for these so- 
lar metallicity progenitors between our prescription and 
this new analysis by Fryer et al. are within ~ 10 — 20% of 
each other. The differences between our prescription and 
the Fryer et al. semi-analytic estimates are well within the 
errors in current detailed stellar evolutionary models and 
collapse simulations. The final black hole mass is primarily 
set by the mass of the star at collapse (our most massive 
remnants are produced in stars that collapse directly to a 
black hole). Differences in the pre-collapse star mass ac- 
count for the bulk of the remnant mass differences between 
stellar models by different groups. These differences are 
primarily driven by different wind prescriptions. 

The most massive remnants are formed through direct 
collapse where the entire (or nearly entire) prcsupcrnova 
star ends up under the event horizon, contributing to the 
black hole mass. In those cases, the explosion energy is too 
low to overcome the gravitational potential of an exploding 
star. For example, a Galactic [Z = Zq) massive star with 
an initial mass of Mzams = 100 Mq is found at presuper- 
nova stage as a WR star with a mass of Msn = 11 — 15 Mq 
and a CO core mass of Afco = 8—11 Mq and this star forms 
a black hole with mass of A/bh = 10 — 13.5 Mq. The lower 
masses in the above example correspond to the Hurley et 
al. winds (§2.1) and the higher masses to the new adopted 
set of winds (§2.2). In both cases the prcsupcrnova star 
was massive enough to form a black hole through direct 
collapse. 

As we will show in the next sections, the change of stel- 
lar wind prescriptions does not affect the compact object 
remnant masses for stars with initial mass Mzams < 30 Mq 
(i.e., stars that form neutron stars, and low mass black 
holes). For the highest mass stars, the new wind pre- 
scription results in more massive presupernova objects 
and heavier black holes. Although the details of super- 
nova calculations are still rather uncertain, relevant studies 
(e.g.. Fryer 1999) indicate that stars with very high initial 
masses Mzams > 100 Mq form black holes through direct 
collapse or at least with significant fall back (i.e., most of 
presupernova mass ends up in the black hole). Therefore, 
the maximum mass of a black hole, which is the main sub- 
ject of this study, will predominantly depend only on (i) 
the employed stellar models and (ii) the adopted set of 
stellar winds, both which set presupernova mass. 

3. RESULTS 

3.1. Standard Prediction 

Here we describe the predictions for the new adopted 
wind mass loss rates and we compare them with the 
previously employed prescription. The results for both 
new and old winds are presented in Figure [1] and in Ta- 
ble [TJ In Figure [T] we present the initial-remnant mass 
relation for three different metallicities; Z = Zq = 0.02. 
Z = 0.3 Zq ^ 0.006 and Z = 0.01 Zq = 0.0002. The 
initial-remnant mass relation shows the compact object 
remnant mass (Mromnant: either neutron star or a black 
hole) for a given initial {M^anis' at Zero Age Main Se- 
quence) star mass. 



For high (Galaxy-like; Z = Zq = 0.02) metalhcity, 
neutron star formation begins at A/zams = 7.7 Mq with 
low mass NSs (A/ns = 1.26 Mq) formed through elec- 
tron capture supernovae (e.g., Podsiadlowski et al. 2004), 
while for higher initial masses (over Afzams = 8.3 Mq; 
see Tab. [T|) NSs form through regular core collapse. In 
a rather wide range 8 <^ A/zams ^ 18 Mq NSs are formed 

with Mns = 1.36 Mq, then for 18 ^ Afzams ^ 20 Mq 

NSs are formed with Afns = 1-86 Mq. There is a bi- 
modal final FeNi core mass distribution due to the mode 
of CO burning - convective or radiative - prior to core- 
collapse; Timmes et al. (1996). The latest results from 
the KEPLER code do not show these same bimodal effect 
on the FcNi core mass, so this effect may not be real 
(Zhang et al. 2008). The detailed study of an alter- 
native neutron star formation mass is underway (Fryer 
& Belczynski, in preparation), we just note that the de- 
tails of neutron star formation do not play crucial role in 
conclusions derived in this study. For Afzams ^ 20 Mq, 
fall back is expected to occur and this rapidly increases 
the mass of the remnant. Depending on the adopted 
limit for maximum NS mass (Afns, max), black hole for- 
mation starts at Afzams ^ 20 Mq (Afns. max = 2.0 Mq) 

or at Afzams - 21 Mq (Af„s,max = 3.0 Mq). Note that 
the NS/BH transition is almost insensitive to the value 
adopted for the maximum NS mass as the initial-remnant 
mass relation rises very steeply for the relevant remnant 
masses (Af remnant > 2 Mq). The steepness of the relation 
is due to the increasing contribution of fall back in the 
final mass of the remnant. 

For lower metallicities, the general features of the initial- 
remnant mass relation for neutron stars are very similar to 
the ones described above, with a natural downward shift of 
initial masses, since lower metallicity stars lose less mass 
in winds. In particular, for very metal-poor environments 
(globular cluster-like; Z = 0.01 Zq ^ 0.0002) NS forma- 
tion starts at Afzams = 6.0 Mq and the transition to BH 
formation is found at Afzams = 18 Mq. Additionally, we 
note that whether we use the old or new wind prescription, 
the initial-remnant mass relation is virtually unaffected for 
neutron stars. We therefore note that the wind mass loss 
rates presented by Hurley et al.(2000) do not differ sig- 
nificantly from the 0/B winds presented by Vink et al. 
(2001) for neutron star progenitors (Afzams ~ 8 — 20 Mq). 
We can derive such a conclusion since in the considered 
mass range, LBV and WR-like winds do not operate. 

For black holes various features of the initial-remnant 
mass relation change significantly with metallicity, and ad- 
ditionally the relation is different for the two sets of winds 
used. For Z = Zq = 0.02 the relation rises steeply for 
Afzams ^ 20 — 24 Mq, and at the high end of this range 
there is a small dip after which the relation resumes its 
rise but at a slower rate. This dip (a WR dip) corre- 
sponds to the onset of helium star winds; the threshold 
mass for removing the entire H-rich envelope is found at 
Mzams. WR = 24.2 Mq and above this mass the stars are 
subject to strong WR winds. At Afzams, lbv ~ 50 Mq we 
observe the onset of very strong LBV mass loss, and as the 
rates differ between old and new winds, the final-initial 
mass relations look different above the LBV formation 
threshold. For the old rates, stronger LBV winds first lead 
to a decrease (a LBV dip) and then a rather slow increase 
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of BH mass with the initial star mass. For the new winds, 
the LBV winds are not as strong, and instead of a LBV dip 
there is only a flattcning-out at A/zams ^ 50 Mq and then 
the BH mass slowly increases with initial mass. The differ- 
ence in relative strength of LBV winds leads to noticeable 
(although not very drastic) differences in black hole mass 
for the two models. Stars with Afzams ^ 100 — 150 Mq 
form Mbh.max ~ 10 Mq BHs in old models (high LBV 
mass loss), while they form Mbh.max ^ 15 Mq BHs in new 
models (moderate LBV mass loss). Note that at solar 
metallicity WR winds are the same for both models. 

For Z = 0.03 Zq = 0.006, BH masses are larger than 
in the case of solar metallicity, as here winds are lower for 
smaller Z so there is a larger mass reservoir at the time of 
BH formation. There are also significant qualitative dif- 
ferences in the shape of the relation. An LBV dip occurs 
for low initial mass (Mzams,LBV ~ 34 Mq), and in partic- 
ular for old winds it is found even before the onset of the 
WR transition/dip (Mzams.WR = 37.4 Mq), while for new 
winds it is found almost at the exact same place as the 
WR dip (Mzams.WR = 33.8 Mq; see Table 1). The shift 
of the onset of WR-like winds corresponds to decreasing 
mass loss with decreasing Z; therefore a higher initial stel- 
lar mass is required to form a naked helium star. Since at 
this Z stars do not lose so much mass, they are more lumi- 
nous and therefore they can reach the LBV phase at lower 
initial masses, so the LBV formation threshold moves the 
opposite way: to lower masses. For high initial masses, 
the BH mass becomes significantly different for the two 
employed models. In particular, for the old set of winds, 
stars with M^ams ~ 25-35 Mq and with M^ams > 100 Mq 
form the most massive BHs: Mbh.max ^ 15 Mq. With the 
new winds, stars with M^ams <; 100 Mq form maximum 
mass BHs: Mbh,max ^ 30 Mq. This increasing difference 
in Afbh,max is due mostly to the dependence of WR winds 
on metallicity. For old winds, WR winds are assumed to be 
independent of Z and therefore the increase in maximum 
BH mass is rather moderate. With the new prescription, 
WR winds are lower for lower Z and we note a rather 
significant increase in maximum BH mass. To give an ex- 
ample, a star with an initial mass of M^ams — 120 Mq 
loses ~ 7 and ^ 25 Mq in WR winds for the new and old 
prescriptions, respectively (more examples in §3.2). 

For Z = 0.01 Zq = 0.0002, the new model shows a new 
(subtle) feature. After the WR dip (at M^ams = 36.9 Mq) 
stars are a subject to strong WR winds, and the BH 
mass increases rather slowly with initial mass. However, 
at Mzams = 92.5 Mq there is a small, but noticeable, 
steepening of the initial-remnant mass relation. At (and 
above) this point stars are sufficiently massive that the 
0/B winds and LBV winds are not strong enough to re- 
move the entire H-rich envelope; so stars never become 
naked helium stars and are not subjected to strong WR 
winds. This is the main reason behind the formation of 
the most massive black holes at such low metallicity. The 
maximum mass BHs are found with Afbh,max 80 Mq 
for Mzams ^ 130 Mq). At higher metallicities, the Z- 
dependent 0/B type winds in conjunction with LBV mass 
loss remove the entire H-rich envelope, and the most mas- 
sive stars are always a subject to WR winds. For the 
old model, the maximum BH mass is found to peak at 
Mbh.max ^ 25 Mq for Mzams ~ 25 — 30 Mq. Iu tliis mass 



range, stars are not subjected to strong WR or LBV mass 
loss, but only to weaker Z-dependent winds (high BH mass 
since the winds are low for low Z; eq. [3]). However, for 
higher initial masses stars lose most of their mass due to 
the effective and Z-independent LBV and WR winds (low 
BH mass since winds are high no matter what Z is used; 
eq. m and eq. [5]) . 

3.2. Examples of Mass Loss 

Here we give some numbers to indicate how much mass 
is lost in specific winds for both models for stars that are 
relevant for maximum black hole mass calculations. 

We start with a star of Mzams = 140 Mq at Z = Zq. 
For Vink et al. winds, this star is Mtms = 63 Mq at the 
end of the main sequence (so it lost Ams = 77 Mq in stel- 
lar winds during its main sequence evolution). When the 
star becomes a WR object its mass is Afwr = 29 Mq (it 
lost Aibv = 34 Mq in LBV winds). At the time of collapse 
the star's mass is M^n = 16.6 Mq (it lost A^r = 12.4 Mq 
in a WR wind) and it forms a (direct collapse) BH of a 
mass Afbh = 15 Mq. For Hurley et al. winds, this star is 
A'^tms = 67 Mq at the end of its main sequence evolution 
(Aijis ~ 73 Mq). When the star becomes a WR object its 
mass is Afwr = 28 Mq (Aibv = 39 Mq). At the time of col- 
lapse the star has Afgn = 11.5 Mq (it lost A^r = 16.5 Mq 
in WR wind) and there is a direct collapse to a BH with 
Mbh = 10.4 Mq. As said before, for solar-like metallicity 
the results are rather similar for both sets of winds, and it 
can be concluded that most of the mass is lost during main 
sequence evolution (^ 53%). The next phase of evolution 
removes less mass (~ 25%) through LBV winds, while fi- 
nally WR winds deplete the star of the relatively smallest 
(but still very significant) amount of mass (~ 12%). 

For comparison, we use the same initial mass of Afzams = 
140 Mq but we shift to models with the very low metal- 
licity Z = O.OIZq. For Vink et al. winds, this star 
is A/tms = 134 Mq at the end of the main sequence 
(Ams = 6 Mq). The LBV phase sets in shortly after the 
star leaves the main sequence and the star explodes dur- 
ing this phase at M^n = 86 Mq (Aibv = 48 Mq). Note 
that the star did not become a WR object (i.e., the H- 
rich envelope was too massive to be removed by prior 
winds). The (direct) collapse leads to the formation of 
a BH with 7\fbh = 78 Mq. For Hurley et al. winds, this 
star is A/tms = 131 Mq at the end of the main sequence 
(Ams = 9 Mq). When the star becomes a WR object its 
mass is Af^r = 58 Mq (Aibv = 73 Mq). Note that in this 
prescription both main sequence winds and LBV winds 
are stronger and the star loses its entire H-rich envelope. 
At the time of collapse, the star's mass is Mgn = 16 Mq 
as it lost an additional Awr = 42 Mq in a WR wind and 
there is a direct collapse to a BH with Mbh = 14.4 Mq. 
The evolution is qualitatively different for both models be- 
cause stronger winds in the Hurley et al. prescription allow 
for the formation of WR stars independent of metallicity, 
while for the Vink et al. winds stars avoid this phase at 
low metallicity and therefore they retain more mass and 
form more massive BHs. It is noted that for such low Z 
only ~ 5% of the star's mass is lost in winds during the 
main sequence, while LBV mass loss is estimated at the 
level of - 34% for Vink et al. and - 52% for Hurley et al. 
Additionally there is a WR mass loss of ~ 30% but only 
for the Hurley et al. prescription. 
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We present here one more model at Z = O.OIZq for 
A^zams = 30 Mq as at this initial mass there is a peak in 
black hole mass (see §3.1) for the Hurley et al. winds. 
This star is Mtms = 29.8 M0 at the end of the main se- 
quence (A„is = 0.2 M0). After the main sequence the star 
never enters LBV nor the WR stage and at the time of col- 
lapse it has Mgn = 25.8 Mq (it lost Ap^^ = 4 Mq in post 
main sequence evolution). After direct collapse a BH is 
formed with Mbh = 23.2 Mq. Only very little mass is lost 
and most of it is lost during post main sequence evolution 
(~ 13%). For Vink et al. winds for the same initial mass, 
the evolution is very similar with almost no mass loss dur- 
ing the main sequence and with only a small amount lost 
during later evolutionary stages (Apms = 0.5 Mq). 

3.3. Effect of LBV Mass Loss 

In Figure [2] (the top panel) we show the initial- remnant 
mass relation for Vink et al. winds at solar metallicity 
{Z = Zq ~ 0.02) for three different levels of LBV mass 
loss. The top line shows the results for LBV winds of 
(dM/dt)LBV2 = 10-^ M0 yr-i (/ibv = 0.1). The maxi- 
mum BH mass in this model is found to be Mbh, max ~ 
80 Mq, very high as the LBV winds are rather weak and 
the most massive stars never lose their entire H-rich en- 
velope (i.e., WR winds do not turn on) and a significant 
fraction of the star's initial mass is retained until the end 
of nuclear evolution. The middle line corresponds to a cal- 
culation with {dM/dt)i,BV2 = lO"'' Mq yr'^ (/ibv = 1-0) 
and the shape of the initial-remnant mass relation is very 
similar to the reference model that we have adopted (see 
Fig. [T] top panel), but with a higher maximum BH mass 
of Mbh. max ^ 20 Mq. Finally, the bottom line denotes 
a model with {dM/dt)i^BV2 = 10"^ Mq yr"! (/ibv = 10) 
and this relation is almost identical to the one that is calcu- 
lated with our old wind prescription (Hurley et al.; § 2.1), 
although here we use a different set of winds (Vink et al.; 
§ 2.2). The maximum BH mass is Mbh, max ~ 10 Mq. 

We choose to set the strength of LBV winds in such 
a way that we can reproduce the most massive known 
BHs in the Galaxy: ~ 15 Mq (e.g., Orosz 2003; Casares 
2007; Ziolkowski 2008). It is found that for /ibv = 1.5, 
which corresponds to a LBV mass loss at the level of 
{dM/dt)i,BV2 = 1.5 X 10"^ Mq yr-\ the maximum BH 
mass reaches 15 Mq in our model with Z ~ Zq (which 
is appropriate for the Galaxy). As a consistency check, 
we test the same model (same LBV wind strength) for 
Z = 0.3 Zq which is the metallicity of the small star- 
forming galaxy ICIO that hosts the most massive known 
BH in a binary: X-1 - 30 Mq (Prestwich et al. 2007; Sil- 
verman & Filippenko 2008). As we can see in Figure[l]the 
model for ICIO can explain BH masses up to ^ 30 Mq. 

The LBV winds were estimated at the level of ~ 10~^ — 
-10-5 Mq yr-i by Vink & de Kotcr (2002). Therefore, 
our choice of {dM / dt)LBY2 = 1-5 x 10"'' Mq yr"^ is found 
very close to the high end of rate range derived by Vink & 
de Koter (2002). However, we demonstrated that for LBV 
mass loss rates as low as <^ 10~^ Mq yr~^, black holes 
would reach unrealistically high masses ^ 80 Mq even for 
high (solar; Z = Zq = 0.02) metallicity (see Fig. O top 
panel). Although, such high mass BHs can not be ex- 
cluded to exist in the Galaxy, we choose to use the masses 
of the known BHs to calibrate our model. If more mas- 



sive BHs are observed then it would be required to lower 
our predicted LBV mass loss rates. We also speculate 
that the discrepancy may originate from the fact that the 
predicted rates do not take into account sporadic bright- 
enings/eruptions of Luminous Blue Variables. If, during 
these episodes, some extra mass is lost then obviously 
the predictions that account only for line/radiation driven 
winds (as presented by Vink & de Koter 2002) might be 
underestimated. For the extreme case of rj Carinae the 
amount of mass lost has been suggested to be substantial 
(^ 10 Mq; Smith et al. 2003). However, it is as yet not 
clear how much mass (if any) is lost in LBV giant eruptions 
in general (e.g., van Genderen 2001). 

We note that in our prescription for LBV mass loss, 
we assume the mass loss does not depend on metallicity. 
However, if most of the mass is lost via metal line-driving, 
we would expect there to be some dependence on metallic- 
ity. It is plausible that such a metallicity dependence for 
the case of LBV winds is weaker than that of OB stars, 
as LBVs are closer to the Eddington limit, which plays 
an important role in the rate of mass loss (see Vink & de 
Koter 2002 for a more detailed discussion). It is clear that 
we need a better understanding of LBV mass loss before 
we can make any definitive conclusions with respect to the 
metallicity dependence of LBV mass loss. 

We also show the results for the Vink et al. winds but 
with slightly lower LBV winds. In Figure [3] (the mid- 
dle panel) we adopt (dM/di)LBV2 = 10~^ Mq yr~^ and 
it is found that the maximum BH mass is: Afbh.max 
100, 35, 20 Mq for metallicities of Z = 0.001, 0.3, i.O Zq, 
respectively. This calculation naturally results in slightly 
higher BH masses, since the LBV winds are less effective 
in removing mass from stars. Note that this model can ac- 
count for the potentially high mass estimate of the Galac- 
tic system Cyg X-1: Mbh = 20 ± 5 Mq, as well as the 
upper mass-range estimates for the black hole in ICIO X-1 
Mbh ~ 34 Mq (e.g., Silverman & Filippenko 2008). 

3.4. Wolf-Rayet winds 

Nugis & Lamers (2000) have presented a study of WR 
wind mass loss rates and their dependence on WR star 
composition. Since, in the stellar models we employ (Hur- 
ley et al. 2000) the actual composition of naked helium 
stars is not available (there is only information on initial 
star metallicity), we can only test the limited formula from 
Nugis & Lamers (2000) that does not depend on the actual 
WR star composition 

log(rfM/(ii)NL = -5.73 + 0.88 log M. (10) 

In Figure [5] (the bottom panel) we show the initial- 
remnant mass relation for our new adopted reference 
model (Vink et al.) in which we have replaced WR star 
winds (cq.[51 Hamann & Koesterke 1998; Vink & de Koter 
2005) with the rates given in eq [TU] adopted from Nugis & 
Lamers (2000). The results are presented for three dif- 
ferent metallicities Z = 1.0, 0.3, O.OIZq. Although now 
WR winds do not depend on metallicity, there is still a 
metallicity dependence on H-rich winds (see eq. [6]and [7]). 
The shape of these initial-remnant mass relations shows 
almost no differences with the new reference model re- 
sults (see Fig. [T|). In fact, the calculations using Nugis & 
Lamers (2000) WR winds resuh in A'/bh.max ^ 30, 80 Mq 
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for low metallicities {Z = 0.3,0.01Zq, respectively) - the 
same values as found in the new adopted model. For high 
metallicity (Z = Zq) the maximum BH mass is found 
to be Afbh.max ^ 20 Mq as compared to ^ 15 Mq for 
the new model. If, in the model with Nugis & Lamers 
(2000) WR winds, we required maximum mass to decrease 
to ^ 15 Mq (maximum mass of known Galactic BHs) then 
we would have to adjust our calibration of LBV winds to 
(dAf/dt)LBV2 = 7 X 10-4 Mq yr-1. 

4. DISCUSSION 
4.1. Comments on Black Hole Masses 

A black hole in the Galactic system GRS 1915 (Mbh ~ 
14 ± 4 Mq) is found to exist in a binary with a low mass, 
evolved stellar companion (Mopt = 1.2±0.2 Mq; K/M III; 
e.g., Greiner et al. 2001) which is filling its Roche lobe. If 
the companion star was born with such a low mass it could 
not possibly have significantly increased the mass of the 
black hole via mass transfer. Belczynski & Bulik (2002) 
proposed an evolutionary scenario in which a companion 
star is born with a low mass and only after the black hole 
is formed it expands during its red giant evolution to fill 
its Roche lobe. In this scenario the observed/current mass 
of the black hole is virtually identical to the mass of the 
black hole at the time of its formation. However, it was 
also argued that the initial mass of the companion may 
have been as high as 6 Mq and therefore, the black hole 
could have increased its mass by ^ 4 Mq (Podsiadlowski, 
Rappaport & Han 2003). 

Another Galactic binary Cyg X-1 is host to a black 
hole (Mbh = 16 ± 5 Mq), and its massive companion 
(Mopt ~ 30 Mq) is estimated to be close to, but not 
quite yet, filling its Roche lobe and is still on the main 
sequence (e.g., Gies & Bolton 1986). The wind mass 
loss rate from the companion is measured at the level of 
2.6 X 10-s Mq yr-i (Gies et al. 2003). If a BH progenitor 
is a very massive star (~ 100 — 150 Mq) then its lifetime is 
^ 4 Myr (e.g. Hurley et al. 2000), and this sets the mini- 
mum age of the system. The main sequence lifetime of the 
companion is estimated to be ^ 5 Myr; we have assumed 
an initial mass of the companion A/opt. zams ^ 40 Mq as it 
has lost about ~ 10 Mq in stellar winds over the course 
of the system's lifetime. This means that the black hole 
may have been accreting from its companion for only 1 
Myr. If we limit the accretion to the Eddington criti- 
cal rate (which is - 10"'^ Mq yr'^ for a 15 Mq BH), the 
black hole in Cyg X-1 may have accreted only 0.1 Mq. 
Had we relaxed the above assumption and allowed this 
BH to accrete at a much higher rate (e.g., Abramowicz et 
al. 1988; King 2002; Ohsuga 2007), the entire mass lost 
in stellar winds by the companion over 1 Myr is 2.6 Mq. 
Only a fraction of this mass can be accreted by the black 
hole, even if the wind is focused in the orbital plane as 
suggested by some observational evidence (e.g., Herrero 
et al. 1995; Miller et al. 2005). This means that the 
black hole mass currently observed in Cyg X-1 is close to 
its formation mass. This must be true if the prior evolu- 
tionary history of Cyg X-1 does not include mass transfer 
via RLOF, and the only mass transfer proceeded through 
stellar winds from the companion (e.g., Ziolkowski 2005). 
Mass transfer via RLOF onto the BH in Cyg X-1 was con- 
sidered by Podsiadlowski et al. (2003). It was concluded 



that if such a phase actually occurred, the black hole mass 
remained basically unchanged due to the fact that RLOF 
was proceeding on a very fast (thermal) timcscale which is 
set by the massive companion, and most of the transferred 
mass was eventually lost from the system. 

The most massive known extragalactic black hole 
(Afbh ^ 30 Mq) is in a binary system ICIO X-1 with a 
massive WR star companion: Mopt ^ 17 — 35 Mq with a 
mass at the high end of this range being the most likely 
case (Silverman & Filippenko 2008). The WR star is well 
within its Roche lobe and the wind mass loss rate was es- 
timated to be at the level of 10"^ - 4 x 10"*^ Mq yr'^ 
(Clark & Crowther 2004). Due to the short lifetime of 
such a massive WR star 0.5 Myr; e.g., Bulik, Belczyn- 
ski & Prestwich 2009), and since only a fraction of wind 
mass can be captured by the black hole, the increase in 
black hole mass due to wind accretion cannot exceed a 
few solar masses. Past evolution of this system has not 
yet been studied in detail and we cannot exclude the pos- 
sibility that there was a phase with RLOF mass transfer 
from the unevolved companion (before it became a WR 
star). Such a mass transfer episode may have, in prin- 
ciple, significantly increased the black hole mass, and a 
thorough investigation of possible initial configurations of 
this system in conjunction with evolutionary calculations 
will be needed to clarify this issue (Valsecchi et al., in 
preparation). 

In conclusion, it appears that the mass we have used 
to calibrate the maximum black hole mass for single star 
evolution in a Galactic-like environment (high metallicity; 
Z = Zq): 15 Mq is consistent with the formation mass 
of Cyg X-1, which hosts one of the most massive BHs 
known in the Galaxy. With such a method, calibrated us- 
ing the most recent metallicity-dependcnt wind mass loss 
rates, we predict that the maximum black hole mass (at- 
tained via single stellar evolution) is Mbh, max ~ 30 Mq 
for an environment of moderate metallicity, like that of 
the ICIO galaxy {Z = 0.3Zq). Our prediction is rather 
striking as the mass of the most massive known BH - that 
which resides in ICIO in binary X-1 - is Mbh = 23 — 34 Mq. 
However, we cannot exclude the possibility that the BH in 
ICIO X-1 has accreted a significant amount of mass from 
its binary companion. 

We also note that if our calibration of LBV winds 
is slightly modified (decreased) to {dM / dt)i^-QW2 = 
10^'^ Mq yr^^, then we predict black holes (single stel- 
lar evolution) to form with masses up to Mbh, max 
20, 35, 100 Mq for Z = Zq (the Galaxy), Z = 0.3 Zq 
(ICIO) and Z = 0.01 Zq (globular clusters), respectively 
(see Fig. [51 middle panel). Therefore, our models can ex- 
plain even the high end of the mass range estimates for 
the most massive known black holes. 

4.2. General Remarks 

In this study, we have analyzed two sets of stellar wind 
mass loss rates that are used in stellar evolution and pop- 
ulation synthesis studies. In particular, we have adopted 
a revised and updated set of stellar winds for the popu- 
lation synthesis code StarTrack. These updates will also 
be included as an option in the BSE population synthesis 
code (Hurley et al. 2002). The previous stellar winds were 
adopted from the compilation of Hurley et al. (2000). The 
new adopted formulae employ metallicity dependent 0/B 
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star winds from Vink ct al. (2001), reduced (clumping) 
Wolf-Rayet winds from Hamann & Koesterke (1998) with 
metallicity dependence from Vink & de Koter (2005) and 
Luminous Blue Variable mass loss rates that were cali- 
brated in this study in such a way that stars could repro- 
duce the most massive known black holes. In fact, our 
primary objective was to check whether within a given set 
of stellar evolutionary models and for available wind mass 
loss predictions, stars can form black holes that are as 
massive as those observed in our Galaxy (~ 14 ±4 Mq for 
GRS 1915 and ~ 16 ± 5 Mq for Cyg X-1) and in external 
galaxies (~ 23 - 34 Mq for ICIO X-1). 

We have demonstrated that for solar metallicity envi- 
ronments (like our Galaxy) , both sets of winds - the pre- 
viously employed (Hurley et al.) and the newly adopted 
(Vink ct al.) - provide remnants massive enough to explain 
known Galactic black holes. However, the predictions are 
quite different for intermediate and low metallicities. In 
summary. Hurley et al. winds remove too much mass from 
stars and within these models one does not expect the for- 
mation of black holes more massive than '^25 Mq even 
at very low metallicities. This is due to the combined ef- 
fects of lack of metallicity dependence for WR winds and 
the rather strong LBV winds adopted in the Hurley et al. 
formulae. Quite contrary, for the Vink et al. winds, maxi- 
mum black hole mass increases with decreasing metallicity, 
and for very low metallicity environments {Z = O.OIZq) 
stars can from black holes as massive as ^ 80 Mq . In par- 
ticular, for a metallicity {Z = 0.3Zq) of the galaxy ICIO 
which hosts the most massive known stellar black hole, our 
newly adopted model predicts the maximum black hole 
mass of Afbh,max = 30 Mq, which is consistent with the 
observed mass. 

Our new results bring the remnant mass predictions of 
StarTrack more in line with the predictions of Fryer & 
Kalogera (2001) and Heger et al. (2003). These studies 
argued for a strong metallicity dependence in agreement 
with the results by many groups studying massive stel- 
lar models: Maeder & Meynct (2001), Meynet & Maeder 
(2003), Heger et al. (2002), Zhang et al. (2008). 

To summarize our results we show the dependence of 
maximum black hole mass on metallicity in Figure [3] for 
both sets of winds. The fact that the maximum black hole 
mass increases with metallicity due to metal-dependent 
winds was recently investigated by stellar evolution com- 
putations with metallicity-dependent main sequence and 
Wolf-Rayet mass loss by Eldridge & Vink (2006) where 
it was found that in lower metallicity galaxies black holes 
are expected to have larger masses. Our results agree and 
to quantify them we provide an approximate formula for 
the maximum black hole mass: 

1403 (^Y - 548^ + 82.9 ^ < 0.11 

-23.7 +75.l(^) - 84.4^ -f 48 ^ > 0.11 

(11) 

where the black hole mass is expressed in [ Mq] and the re- 
lation holds true for our newly adopted set of stellar winds 
(Vink et al.). These formulae are only valid for metallic- 
ities 0.01 < (Z/Zq) < 1.5; the range is set by the limits 



of the employed stellar evolutionary models (Hurley et al. 
2000). The maximum black hole mass rises from ~ 10 Mq 
for super-solar metallicity (Z — 1.5Zq) to ~ 80 Mq for 
extremely low metallicity (Z — O.OIZq). The Eddington 
limited accretion onto an 80 Mq black hole may give a rise 
to high X-ray luminosity ^ 10^° erg s"""^ and ~ 2 x 10"*° 
erg in case of hydrogen and helium accretion, respec- 
tively. Such high luminosities are comparable with X-ray 
emission of many known ULXs (e.g., Madhusudhan et al. 
2008; Zampieri & Roberts 2009). 

For the most massive stars that will form maximum 
mass BHs, we have estimated the contribution of specific 
winds in the overall mass loss from a star. For high metal- 
licities, most mass is lost during main sequence evolution 
(~ 50%), later LBV winds also remove a significant part 
of the star's initial mass (~ 25%) while finally the exposed 
naked helium star further loses mass (~ 10%) in WR-like 
winds. For very low metallicities, the majority of mass is 
lost during the LBV phase (^ 30%) with only small mass 
loss on the main sequence (^ 5%), while the WR phase 
is not encountered (i.e., star never loses its H-rich enve- 
lope). For intermediate metallicities, the mass is lost in 
comparable quantities during the main sequence and LBV 
phases (^ 30%), with a small contribution to mass loss 
from WR-hke winds 5%). 

Our approach, based on one set of stellar evolution mod- 
els by Hurley et al. (2000), matches observational con- 
straints on black hole masses. However, stellar evolution 
models continue to produce very different results. The 
differences in evolution models produce uncertainties in 
stellar mass, core mass and thus in final black hole mass. 
To illustrate these differences, we have compared our re- 
sults based on the Hurley et al. (2000) formulae to models 
obtained with the EZ stellar evolutionary code (Paxton 
2004) which we modified to include our specific wind mass 
loss prescriptions. The results of our comparisons are pre- 
sented in Table [2] and [31 We note that the results, both 
for the analytic formulae and for detailed evolutionary cal- 
culations, show the same trends; for example the mass of 
the He core at the end of Main Sequence is higher for the 
old Hurley et al. winds as compared to models with the 
new Vink et al. winds. But the EZ code typically pre- 
dicts He core masses that are ^ 20% higher than those 
produced by our analytic formulae. Our analytic formulae 
were based on the revised Eggleton evolutionary code (Pols 
et al. 1998), while the EZ code is also based on the original 
Eggleton code but revised by Paxton (2004). Even with 
very similar codes, properties in the core still have large 
differences. The shape of our black hole mass distribution 
is sensitive to these uncertainties. However, the maximum 
black hole mass is produced in systems where the entire 
star collapses down to a black hole. To first order, the max- 
imum black hole mass depends upon the final stellar mass 
and not on the internal stellar structure. The final stellar 
mass depends mostly on the wind mass loss rates, but it 
changes also with the adopted stellar evolutionary model 
(which results in different luminosities, radii, etc. so the 
same mass loss prescription gives a different rate). Thus, 
our results are in the end limited by the details of stellar 
evolutionary modeling uncertainties. It would be very use- 
ful if the stellar evolutionary community could provide the 
final stellar masses that are obtained with different codes. 
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Such a comparison would set the uncertainty of our resuh 
on the final black hole mass. One can hope that the final 
stellar masses obtained with various detailed evolutionary 
codes but with the same set of wind mass loss rates would 
not differ by more than ~ 30%. 

Some stellar evolutionary models employ rotation 
(Heger et al. 2000; Maeder & Meynet 2001; Meynet & 
Maeder 2003; Vazques et al. 2007; de Mink et al. 2009). 
For stars that rotate fast (either low metallicity single stars 
or tidally locked stars in close binaries), the effect of ro- 
tation may play an important role on their evolution and 
subsequent formation of the compact object remnant. The 
primary effect of rotation is mixing of additional material 
into the central part of the star where burning is taking 
place, thus leading to the formation of a more massive 
core (as compared with non- rotating models). This should 
lead, in principle, to the formation of more massive rem- 
nants (e.g., it should increase the maximum black hole 
mass). However, rotation makes a star oblate, increasing 
the temperature of the polar regions. Thus, there is more 
radiation to drive winds, and the net effect may be an in- 
crease in mass loss from a rapidly spinning star leading 
to a lower presupernova mass and so, reduction of the fi- 
nal compact object mass. It is not completely clear what 
would be the full effect of rotation on the maximum black 
hole mass, especially since the initial rotation of massive 
stars is not known, and this is particularly true for low 
metallicities. 

Finally, we want to emphasize that our conclusions arc 
obtained only for single star models. In other words we 
do not take into account potential effects of binary evolu- 
tion. Those, in principle, may increase the mass of a black 



hole (via accretion from a close companion) or reduce its 
mass (by mass loss from its progenitor). Additionally, the 
black hole masses that are known come from the obser- 
vations of very specific types of sources: X-ray binaries; 
and thus it may be possible that either some evolutionary 
and/or observational selection effects hide the true, intrin- 
sic maximum black hole mass. Therefore, by its nature, 
this is only a proof-of-principle study demonstrating that 
with some (widely used) evolutionary models and a set 
of (updated) metallicity-dependent stellar winds, regular 
single star evolution can naturally explain the metallicity- 
dependence of the most massive black holes that are known 
today. Additionally, it needs to be highlighted that all 
three major components used in our study (i) stellar evo- 
lution models, (ii) wind mass loss rate predictions and 
(in) core collapse/supernova calculations that are all very 
important in remnant mass estimates are burdened with 
significant uncertainties and sometimes even unknowns. 
We have attempted to collect some of the most recent and 
widely used results (both observational and theoretical) to 
provide a self-consistent physical model for remnant mass 
calculations. However, we note that this model is still sub- 
ject to future adjustments and changes as more constraints 
will become available. 
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Table 1 

Characteristic Properties of Initial-Remnant Mass Relations' 











Mzams[ Mo]: 








-^remnant • 






WD/NS*^ 


NS/BH 


NS/BH 


NS/BH 


WR 


LBV 


-^^bh.max 


Model 


Z 




(2.0 Mo) 


(2.5 Mo) 


(3.0 Mo) 






[Mo] 


Hurley 


Zq 


7.7 (8.3) 


20.2 


20.8 


21.3 


24.2 


48.5 


10.5 


Hurley 


0.3 Zq 


7.0 (7.7) 


19.2 


19.4 


19.6 


37.4 


34.3 


16.0 


Hurley 


0.01 Zq 


6.0 (6.8) 


18.0 


18.2 


18.3 


36.2 


32.0 


24.2 


Vink 


Zq 


7.7 (8.3) 


20.2 


20.8 


21.2 


24.2 


49.7 


15.0 


Vink 


0.3 Zq 


7.0 (7.7) 


19.0 


19.2 


19.4 


33.8 


33.9 


28.3 


Vink 


0.01 Zq 


6.0 (6.8) 


17.9 


18.1 


18.2 


36.9^^ 


32.0 


79.1 



'^Initial star mass (A/zams) is given for: transition of WD to NS formation and NS to BH 
formation for three assumed maximum NS masses A/ns,max = 2.0, 2.5, 3.0 Mq- Minimum initial 
mass over which Wolf-Rayet star forms and over which Luminous Blue Variable forms. In the 
last column we list the maximum BH mass that is formed in a given model. 

^Number in parenthesis: initial star mass for transition from electron capture supernova NS 
to core collapse NS. 

■^WR stars form only in the limited range of M^ams = 36.9 — 92.5 Mo; for details see§3.1. 



Table 2 

Comparison of Models: M^ams = 50 Mq 



Z /Wind 






' Ills [-^'^-'-j 1 J 


-'^-'bilL ^*-^0J 


0.02/NoWind 


50.0 (50.0) 


22.5 (18.9) 


3.8 (4.2) 


(45.0) 


0.02/Hurley 


42.3 (41.3) 


21.2 (14.6) 


3.9 (4.3) 


(7.9) 


0.02/Vink 


42.3 (40.6) 


20.6 (14.3) 


3.9 (4.4) 


(9.4) 


0.004/NoWind 


50.0 (50.0) 


22.0 (18.9) 


4.2 (4.5) 


(45.0) 


0.004/Hurley 


46.9 (46.5) 


21.6 (17.1) 


4.2 (4.6) 


(6.7) 


0.004/Vink 


47.4 (47.3) 


21.6 (17.5) 


4.2 (4.6) 


(14.8) 


0.0003/NoWmd 


50.0 (50.0) 


22.0 (18.9) 


4.3 (4.6) 


(45.0) 


0.0003/Hurley 


49.2 (49.1) 


22.0 (18.5) 


4.3 (4.6) 


(8.6) 


0.0003/Vink 


49.7 (49.6) 


22.0 (18.7) 


4.3 (4.6) 


(17.9) 



^ Values for detailed evolutionary calculations (and the analytic formulae) are given for a star 
with initial mass Afzams = 50 Mq. Mtms^star mass at the end of Main sequence; Afhc~helium 
core mass at the end of MS; Tms-MS lifetime; Mbh-final mass of a remnant (this is given only for 
our calculations as evolutionary models were not evolved beyond main sequence. For details see 
Sec.4.2. 



Table 3 

Comparison of Models: Mzams = 100 M0 



Z/Wind 




Mhe[ Mq] 


Ti„s[Myr] 


Mbh[ Mq] 


0.02/NoWind 


100.0 (100) 


51.9 (48.3) 


2.7 (3.4) 


(90.0) 


0.02/Hurley 


63.0 (60.3) 


43.9 (24.4) 


2.7 (3.5) 


(9.8) 


0.02/Vink 


60.9 (59.1) 


43.1 (23.7) 


2.8 (3.5) 


(13.7) 


0.004/NoWind 


100.0 (100) 


50.6 (48.3) 


2.9 (3.5) 


(90.0) 


0.004/Hurlcy 


86.2 (80.5) 


47.3 (36.0) 


2.9 (3.5) 


(11.6) 


0.004/Vink 


94.2 (84.0) 


49.3 (38.1) 


2.9 (3.5) 


(33.5) 


0.0003/NoWind 


100.0 (100) 


50.6 (43.4) 


3.0 (3.5) 


(90.0) 


0.0003/Hurley 


96.5 (95.1) 


49.4 (41.0) 


3.0 (3.5) 


(12.4) 


0.0003/Vink 


99.5 (98.2) 


50.5 (42.5) 


3.0 (3.5) 


(42.9) 



Same as Table [5] but for initial star with mass M^ams — 100 Mq. 
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Fig. 1. — The initial- remnant mass relation for single stellar evolution for two wind prescriptions: the previously 
used Hurley et al. winds (§2.1) and newly adopted modified Vink et al. winds (imposed LBV winds at the level of 
1.5 X 10""' M0 yr~^; §2.2). Top panel: Results for solar metallicity that correspond to the stellar field populations in 
Galaxy. The predicted maximum black hole mass Mbh,max 15 Mq for new and ^ 10 Mq for old winds is consistent with 
the most massive stellar black holes observed in our Galaxy (e.g., in GRS 1915 A/bh = 14±4 M0). Middle panel: Results 
for moderate metallicity that correspond to stellar populations in galaxy ICIO which hosts the most massive known stellar 
black hole (Mbh = 23 — 34 Mq). Note that the predicted maximum black hole mass Mbh,max ~ 30 Mq for new winds is 
consistent with the measurement in ICIO, while Mbh. max ~ 15 Mq obtained for old winds appears to be significantly too 
small. Bottom panel: Results for very low metallicity that correspond to stellar populations of Galactic globular clusters 
or metal-poor galaxies. The maximum black hole mass may reach Mbh, max ^ 80 Mq or 25 Mq for the new and old 
wind prescriptions, respectively. Note the change of vertical scale from panel to panel. 
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Fig. 2. — The initial-remnant mass relations for alternative parameter choices. Top panel: The relation for Vink et 
al. winds with the Galaxy-like metallicity {Z = Zq) with three different assumed levels of LBV winds: 10~'^ Mq yr~^ 
(/ibv = 0.1, top line), 10"'' Mq yr'^ (/ibv = 1-0, middle line), 10"^ Mq yi"^ (/ibv = 10, bottom line). The maximum 
BH mass is: A/bh.max ^ 80, 20, 10 Mq for the top, middle and bottom lines, respectively. Note that Galactic black 
holes reach masses of ^ 15 Mq, implying the LBV winds at the level of ^ 10~^ — 10~^ Mq yr~^. Middle panel: The 
relation for the Vink et al. winds but with lower (than adopted in the reference model) LBV winds: 10~^ Mq yr~^. 
Note the high maximum BH masses for this model Mbh.max 100, 35, 20 Mg for Z = 0.001, 0.3, I.OZq, respectively. 
Bottom panel: The relation for Vink et al. winds with modified WR winds that were adopted from Nugis & Lamers for 
these calculations. Note that results for low metallicity {Z = 0.01,0.3Zq, upper curves) are almost identical to our new 
reference model, while the high metallicity model results in a higher maximum BH mass (Mbh.max ~ 20 Mq) than the 
reference calculation 15 Mq). 
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Fig. 3. — The dependence of maximum black hole mass on metallicity for our previously used set of stellar winds (Hurley 
et al.; §2.1) and for new adopted winds (Vink et al.; §2.2). Note that the maximum BH mass is similar for solar-like 
metallicities, while it is significantly larger for the new winds in metal-poor environments. 



